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1 Introduction 

It is well known that the characteristic Lie algebra introduced by A.B. Shabat in 1980, plays 
the crucial role in studying the hyperbolic type partial differential equations. For example, if 
the characteristic algebra of the equation is of finite dimension, then the equation is solved in 
quadratures, if the algebra is of finite growth then the equation is integrated by the inverse 
scattering method. More details and references can be found in pQ. Recently it has been 
observed by A.V. Zhiber that the characteristic algebra provides an effective tool for classifying 
the nonlinear hyperbolic equations. Years ago the characteristic algebra has been used to classify 
integrable systems of a special type j2j. However, the characteristic algebras have not yet been 
used to study the discrete equations, despite the fact that the discrete equations had become 
very popular the last decade (see, for instance, survey 

In this paper we show that the characteristic algebra can be defined for any discrete equation 
of the hyperbolic type and it inherits most of the important properties of its continuous counter- 
part. However, it has essentially more complicated structure. The work was stimulated by [I], 
where the discrete field theory is studied and the question is posed which types of algebraic 
structures are associated with the finite field discrete 3D Toda chains. 



2 Invariants and vector fields 

Consider a discrete nonlinear equation of the form 

t{u + 1, v + 1) = f(t(u, v),t{u + 1, v), t(u, v + 1)), (1) 

where t = t(u, v) is an unknown function depending on the integers u, v, and / is a smooth 
function of all three arguments. The following notations are used to shorten formulae: t u — 
t(u + 1, v), t v = t(u, v + 1), and t uv = t(u + 1, v + 1). By using these notations one can rewrite 
the equation Q as follows t uv = f(t,t u ,t v ). Actually, the equation is a discrete analog of the 
partial differential equations. Particularly, the class of equations contains difference schemes 
for the hyperbolic type PDEs on a quadrilateral grid. 
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The notations above are commonly accepted, but not very convenient to indicate the iterated 
shifts. Below we use also different ones. Introduce the shift operators D and D, which act as 
follows Df(u,v) = f(u + l,v) and Df(u,v) = f(u,v + 1). For the iterated shifts we introduce 
the notations fj = D^(f) and fj = Di(f), so that t(u+l, v) =t±, t(u,v + l) = i±, t(u+2, v) = t 2 , 
t(u, v + 2) = t,2 and so on. 

The equation (^Q) is supposed to be hyperbolic. It means that it can be rewritten in any of 
the forms: t u = g(t,t v ,t uv ), t v = r(t,t uv ,t u ), and t = s(t u ,t uv ,t v ) with some smooth functions 
g, r, and s. 

A function F = F(v,t, ti,ti, . . . ), depending on v and a finite number of the dynamical 
variables is called u-invariant, if it is a stationary "point" of the shift with respect to v so that 
(see also j3]) 

DF(v,t,h,t 1 ,...) = F(v,t,ti,ti, ...), (2) 

and really the function F solves the equation F(v + 1, t v , f, /i, £2, ■ ■ ■ ) = F(y, t, t±, ii, . . . ). Exa- 
mining carefully the last equation one can find that: 

Lemma 1. The v-invariant does not depend on the variables in the set {ij}J^ =1 . 

If any v-invariant is found, then each solution of the equation (JTJ can be represented as 
a solution of the following ordinary discrete equation F(t,t u , . . . ,tj) = c(u), where c(u) is an 
arbitrary function of u. 

Due to the Lemma ^ the equation (j2j) can be rewritten as 

F(v + l,t v ,f,f 1 ,...)=F(v,t,t 1 ,t 2 ,...). 

The left hand side of the equation contains t v , while the right hand side does not. Hence the 
total derivative of DF with respect to t v vanishes. In other words, the operator X\ = D^-jr-D 
annihilates the u-invariant F: X\F = 0. In a similar way one can check that any operator of the 
form Xj = D~i -^jr-Di , where j > 1, satisfies the equation XjF = 0. Really, the right hand side 
of the equation D J F(v, t, t\, i±, . . . ) = F(v,t,ti,ti, . . .) (which immediately follows from ©) 
does not depend on t v and it implies the equation XjF = 0. As a result, one gets an infinite set 
of equations for the function F. For each j the operator Xj is a vector field of the form 

Consider now the Lie algebra L v of the vector fields generated by the operators Xj with the 
usual commutator of the vector fields = X^Xj — XjX^. We refer to this algebra as 

characteristic algebra of the equation 

Remark 1. Note that above we started to consider the discrete equation of the form Q con- 
jecturing that it admits nontrivial v-invariant. The definition of the algebra was motivated by 
the invariant. However the characteristic Lie algebra is still correctly defined for any equation 
of the form Q even if it does not admit any invariant. 

3 Algebraic criterion of existence of the invariants 

Theorem 1. The equation (^Q) admits a nontrivial v-invariant if and only if the algebra L v is 
of finite dimension. 
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Proof. Suppose that the equation admits a non-constant solution. Then the following 
system of equations 

XjF(v,t,t u t 2 ,...,t N ) = 0, j>l (4) 

has a non-constant solution. It is possible only if the linear envelope of the vector-fields {Xj}J? =1 
is of finite dimension. 

It is worth mentioning an appropriate property of the vector fields above. If for a fixed j the 
operator Xj is linearly expressed through the operators X\, X2, • • • , Xj—±, then any operator Xk 
is a linear combination of these operators. Really, we are given the expression Xj = a\X\ + • ■ ■ + 

dj-xXj-t. Note that X j+1 = D^XjD, hence X j+1 = D~ 1 (a 1 )X 1 H h D -1 (oj_i)Xj. Thus, 

in this case the characteristic algebra is generated by the first j — 1 operators which are linearly 
independent. Due to the classical Jacoby theorem the system Q) has a non-constant solution 
only if dimension of the Lie algebra generated by the vector fields Xi is no greater than N. 
Thus, one part of the theorem is proved. 

Suppose now that the dimension of the characteristic algebra is finite and equals, say, N, 
show that in this case the equation (pQ) admits a ^-invariant. Evidently, there exists a function 
G(t,t±, . . . , tjv)) which is not a constant and that XG = for any X in L v . Such a function 
is not unique, but any other solution is expressed as h(G). Due to the construction the map 
X — > D~ 1 XD leaves the algebra unchanged, hence G\ = DG is also a solution of the same 
system XG = and therefore G\ = h(G). In other words, one gets a discrete first order 
equation: DG = h(G), write its general solution in the following form: C = F(v,G) where C 
does not depend on v. Evidently the function F found is just the w-invariant needed. ■ 



4 Computation of the characteristic algebra 

In this section the explicit forms of the operators {Xj} are given. We show that the operators 
are vector fields and give a convenient way to compute the coefficients of the expansion (j3J). 

Start with the operator X\. Directly by definition one gets X\F(t, t\, . . .) = D^ 1 -^-F(t v , f, 
fx, . . .). Computing the derivative by the chain rule one obtains 

™--) = fl - , (^ + ^ + ^ + -) F( '" / ' A -- ) ' 



and finally 



So the operator X\ can be represented as X\ — ^i~§t~ where the coefficients Xj are found 

(=0 

as Xj = D^ 1 { ^Q t ~ l ^ for j > and xq = 1. Actually the coefficients can be computed by the 
following more convenient formula 

x j+1 = xtDix^D^xt) ■ ■ ■ Di( Xl ), (6) 

or Xj + \ = XjD 3 {x{). Really, 

j+1 dt v dfj-idfj-2 dt v \dt v J \dt v J dt v 

= D(D j (x 1 )---D(x 1 )x 1 ). 
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To find X 2 , use the following formula 

di + Xl df + X2 df 2 + " ' ) F{tvi f ' /l ' ' ' ' ] - 

After opening the parentheses and some transformation the right hand side of the last formula 
gets the form 

X 2 F(t, t 1 ,...) = D- 1 (x x {f)%- + X^h) A + Xt(f 2 ) A + ■■■) F(t v , 
So the operator can be written as 

X 2 = Z)" 1 (X^/)) A + 5- 1 A + ( Xl (/ 2 )) A + . . . . (7) 
Continuing this way, one gets 

X 3 = A + 5-1 A + 5-1 (^(/a)) A + • • • . (8) 

One can derive an alternative way to compute the coefficients of the vector fields X k above. 
Represent the operators as follows 

00 d 

Xk = ^n kj — . (9) 
3=0 1 

We will show that the coefficients njy of the operators satisfy the following linear equation 

Dn k+W = jy (A) n kJ + D* (A) n kJ+1 + (A) Dn k+1 „ (10) 

closely connected with the direct linearization of the initial nonlinear equation ([T]). In order to 
derive this formula, apply the operator X k to the iterated shift fj = f(tj,tj + i, fj-i) and use the 
chain rule 

XkLfj) = & {%) D-\X k ^{f 3 ^)) + Lfi (A) D' l X k ^{f 3 ) + Dif^Xkifj-MU) 

Comparison of the two representations (JSJ) and @ of the operator X k yields X k (fj) =Dn k+ \j + \. 
By replacing in 1)11(1 X with n one gets the formula 1)10)1 required. 

The characteristic algebra is invariant under the map X — » D~ l XD. First prove the formula 

D~ l X\D = D (x{)X\. To this end use the coordinate representation of the operator X\ = 

00 

^2 Xi-M- and the formula © for the coefficients. Then check that D~ 1 X 2 D = pX 2 . Really, 

i=o i 

D~ X X 2 B = D~ 1 D~ 1 X 1 DD = D~ 1 D~ 1 X 1 DD = D~ 1 (D~ 1 (x 1 ))X 2 . 

Obviously similar representations can be derived for any generator of the characteristic algebra. 
The following statement turns out to be very useful for studying the characteristic algebra. 

Lemma 2. Suppose that the vector field 

00 „ 

3=0 J 
satisfies the equation 

D(X) = X, (12) 
then X = 0. 
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Proof. It follows from ()12|) that 



OO q OO q 

^2D(x j (ii,i-i,t,t 1 ,t 2 , . . .,t k] ))-^— = J ^2x j (ii,i-i,t,t 1 ,t 2 ,. . -,t kj ) — . 



3=0 



+ 1 



3=0 



Comparison of the coefficients before the operators in both sides of this equation yields: 
xq = 0, Xi = 0, . . . and so on. ■ 

5 The commutativity property of the algebra 

One of the unexpected properties of the characteristic Lie algebra is the commutativity of the 
operators Xj. Consider first an auxiliary statement. 

Lemma 3. The coefficients x k i of the vector fields X k , k > 1, i > do not depend on the 
variable t v . 

In other words, the coefficients of the expansions 

°° d 

3=0 3 

satisfy the equation -^-x ki = 0. 

Proof. For the operator X\ one has xij+i = Z)" 1 (J^^J- But the function fj does not depend 

on t 2 = Dt v and on £3, £4, . . . as well. Hence the coefficients do not depend on t\ = t v . Similarly, 
the functions X\{fj) do not depend on t 2 , so x 2 ,j+i = -D -1 (Xi(fj)) do not contain t v . One can 
complete the proof by using induction with respect to k. ■ 

Now return to the main statement of the section. 

Theorem 2. For any positive integers i, j the equation holds [Xi,Xj] = 0. 
Proof. Remind that Xj = D~ i-^-D 3 , so that one can deduce 



D^^-D^D^-^-D 3 



[Xi,Xj] 

Suppose for the definiteness that k = j — i > 1. Then 
[Xi,Xj]=D 

But due to the Lemma |3 the last expression vanishes that proves the theorem. 







D i = D~ { 


— , X k 




dty 



OO , rj 

& = D~ l Y —(x ki )^-D\ 

3=0 j 



6 Characteristic algebra for the discrete Liouville equation 



The well known Liouville equation P ~ = e v admits a discrete analog of the form [5] 



dxdy 



(13) 



6 



I.T. Habibullin 



which can evidently be rewritten as 
1 

Specify the coefficients of the expansions ©-Q) representing the vector fields X% and X2 for the 
discrete Liouville equation (|13|). Find the coefficient x\ = D~ l remind that f(t,t u ,t v ) = 

j(t u - l)(t v - 1) and g(t- v ,t,t u ) = ^tut-v + 1, 

x 1 = D^(?L)=D-^ tu - 1 ^- 9 - 1 - tu 



dt v ) \ t ) t- v t — 1 

Similarly 

x 2 = D 



df dt v J h-lt-l 

J-i t 

It can easily be proved by induction that Xj = Yi j^zi ^ or 3 — ^' here to := t- Remind also that 

k=o k 

xq = 1. So that the vector field is 

x _ d h d t x t 2 d tit 2 t 3 d_ 

1 ~ dt t-ldh (t-l)(ti-l)a* 2 (t- l)(ii - l)(t 2 - 1) 9t 3 '"' 

It is a more difficult problem to find the operator X2, it can be proved that 
X2 = — 1 (-^1 — X— 1), 

[t — V 1)^ — 1) 

where the operator X_i is defined as follows X_i = D-J^D" 1 . For the coefficients x_ 7 - of the 



00 a j ~ l t -1 

expansion X-i = Yl X -JW~ one can deduce the formula = n — , so that the operator 

j=o ' J ' k=o k 

is represented as 

X _ = & I h ~ l 8 I (^-^fe-l) g , (*1-1)(*2-1)(*3-1) g 

1 dt t 9ti tti dt 2 Uit 2 dt 3 

X 3 = D~ l X 2 D = { ^ 1 ~ 1) ^ 1 X 2 . 

(t-2-l)t- 2 

Theorem 3. Dimension of the characteristic Lie algebra of the Liouville equation (|13[) equals 
two. 

Proof. It is more easy to deal with the operators Y + = (t — 1)X\ and Y_ = tX_i rather than 
the operators X\ and X-\. In order to prove the theorem it is enough to check the formula 

[Y+,y_]-y+ + Y_ = o. (14) 

Denote through X the left hand side of the equation (|14jl and compute D(X) to apply the 
Lemma |2 It is shown straightforwardly that 
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Compute now the shifted operator D([Y + , Y_]) = [D(Y + ), D(Y_)] which can be represented as 
follows 



[D(Y + ),D(Y_ 



ti-1 

h 



dtj'h-i 



-4 



A x + A 2 + A 3 + A A , 



where 



A, 



U - 1 



A, 



(*!_!)(* -1)0 tj 



-y_ 



tt a a 

J4 



tl ^'tl-iat 

Direct computations give 

A! = [r + ,y_]-^y_ 



U-l UY, 



A a 



An 



ti dt h-1 
(ti - l)(t - 1) d tii 5 



1 - t - ti 
^ dt : 



dt' h-ldt 



A, 



1 - t - ti 5 
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Al = — 



ti-1 dt' "™ 1 dt' 
Summarizing all computations above one gets the result 



D[Y + ,Y_] = [Y + ,Y_] 



1 



ti-1 



-y_ 



i 

/a 



y+ + i + 



t - 1 



/i 



+ 



a 



h-lj dt' 



which implies D([Y + , V_] — Y" + +Y1) = [y + , Y_] — y + +y_. Now apply Lemma El to the function 

x = [y+, y_] - y+ + y_ to get z = o. ■ 



7 How to find the invariants? 

In this section two examples of the equations with the invariants are shown. Start with a simple 
one. 

Example 1. Consider a linear equation of the form 

tuv — t u -\- t v t -\- 1 

so that f(t, t u ,t v ) = t u + t v — t + 1 and g(t_ v ,t, t u ) =t u — t + t- v — 1. It is easy to see that the 
algebra is of one dimension 

d d d X — 

1 dt dt u dt2 2 

The first integral of the least order for the equation X\F = can be taken as Go = ti — t. 
Evidently it solves the equation DGq = Gq + 1. Now, the invariant is to be taken as I = t\—t — v, 
because DI = t uv — t v — v — 1 = t u — t — v = I. Any other ^-invariant of the equation can be 
represented as F = F(I, DI, D 2 I, D k I). 

Example 2. Return to the discrete Liouville equation discussed above. Find the intersection 
of the kernels of the operators X\ and X-\, hence they constitute the basis of the characteristic 
algebra of the Liouville equation (|13j) . To this end first solve the equation X\F = which is 
reduced to the following infinite system of the ordinary differential equations: 

dt dti dt2 



1 ti/(t-l) tit 2 /(t-l)(ti-l) 
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So that the invariants of the vector field X\ are Iq = I\ = D(Iq), I 2 = D(I{), .... Change 
the variables in the vector fields by setting t = t, t% = Iq, t 2 = lit ■ ■ ■ ■ Then one gets 

v d d h + i a t 2 + i d 

X\ = — and A_i = -7— 



dt dt t(t - 1) dh ht(t - 1) dt 2 

Now solve the equation & = jjvjzl— to find the common solution Fo = (r^zj + lj (^j^- + 1 
of both equations X\F = and X_\F = 0. The invariant is 

/(*,*!, t 2 ) = F (t,h,t 2 ) = (-^— + 1 ) ( ^ + 1 



n *i — 1 y V * 

(see also Evidently, the basis of common invariants can be chosen as follows: Fq,D(Fq), 
D 2 (Fq), .... It can be easily verified by the direct computation that Fq solves the equation 
DF = F Q . 

Corollary 1. In the case of the Liouville equation each solution of the system of the equations 
X\F = ; X-±F = gives the v-invariant. 

Proof. Each solution of this system can be represented as F = F(Fq, D(Fq), D 2 (F), . . . ), hence 

DF(F , D(F ),D 2 (F), . . . ) = F(D(F ), D(D(F )), D(D 2 (F)), . . . ) 

= F(F ,D(F ),D 2 (F),...). U 

8 Discrete potentiated Korteweg— de Vries equation 

Discrete equation of the form 

(t uv - t + p + q){t u - t v + q - p) = q 2 - p 2 (15) 

is called the discrete potentiated KdV equation (DPKdV). The constant parameters p, q are 
removed by a combination of the shift and scaling transformations of the form t(u, v) = ht(u, v) — 
pv—qu, where h 2 = q 2 —p 2 \iq 2 —p 2 > and h 2 = —q 2 +p 2 \iq 2 —p 2 < 0. After this transformation 
the equation (fT3|) takes one of the forms 

The algebraic properties of the equation do not depend on the choice of the sign before the unity, 
we will take the upper sign. Study the problem how to describe the characteristic Lie algebra 
for the equation 

(tuv-t)(t u -t v ) = 1. (16) 

Represent the equation ((16]) as t uv = f(t,t u ,t v ), where f(t,t u ,t v ) = t+ t J_ tv ■ Set x = = 
77 — h-^i, then x = (f — t) 2 . Below we use the upper index to denote the shift with respect to v 

and lower index to denote the shift with respect to u so that D~ k Di x = X:. Then due to the 
general rule the operator X\ is written as 

X 1 = — + x 1 — + x^l— + ■■■ + x l x\ ■ ■ ■ x)— + • • • . 

at ot\ ot2 J otj+i 

Theorem 4. Characteristic Lie algebra of the DPKdV equation (fTB|) is of infinite dimension. 
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Proof. In the case of the DPKdV equation one gets = 1, = —x, = x. The for- 
mula (|1U|1 is specified 

Dn k+ i J+1 = n k j - Xjn k j + i + XjDn k+ ij. 

It is more convenient to write it as 

n k+1J+1 = D^rikj - XjD^rikj+x + x]n k+1 j, 

remind that D~ 1 xj = xj. It is clear now that the coefficients are polynomials of the finite 
number of the dynamical variables in the set {x* }^ =0 . Order the variables in this set according 
to the following rule: ord(x") > ord(Xp) if j > p and ord(xy i ) > ord(Xp) if m > n. 

Lemma 4. For any positive k,j the coefficient n k j for the operator X k in the DPKdV case is 
represented as 

k 

n k,j = ni'k,jXj_± + r k j 

where m k j and r k j are polynomials of the variables with the orders less than the order of Xj_ lr 
moreover, m k j is a monomial. 

Lemma can be easily proved by using the induction method. It allows to complete the proof 
of the theorem. The vector fields {X k } k *L 1 are all linearly independent. ■ 

Due to the theorem the DPKdV equation do not admit any u-invariant. It is not surprising, 
because the equation can be integrated by the inverse scattering method or, in other words, it is 
S-integrable, it is well known in the case of partial differential equations that only C-integrablc 
equations (Liouville type) admit such kind objects called x- and y-integrals. 

9 Conclusion 

The notion of the characteristic algebra for discrete equations is introduced. It is proved that the 
equation is Darboux integrable if and only if its characteristic algebras in both directions are of 
finite dimension. The notion can evidently be generalized to the systems of discrete hyperbolic 
equations. It would be useful to compute the algebras for the periodically closed discrete Toda 
equation, or for the finite field discrete Toda equations found in 0] and [S], corresponding in the 
continuum limit to the simple Lie algebras of the classical series A and C. 
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